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ABSTRACT 

In the relativistic plasma surrounding a pulsar, a subluminal ordinary-mode 
electromagnetic wave will propagate along a magnetic field line. After some 
distance, it can break free of the field line and escape the magnetosphere to 
reach an observer. We describe a simple model of pulsar radio emission based 
on this scenario and find that applying this model to the case of the Vela pulsar 
reproduces qualitative characteristics of the observed Vela pulse profile. 

Subject headings: pulsars: general — pulsars: individual (Vela) 

1. Introduction 

Soon after the discovery of pulsars, Goldreich & Julian (1969) deduced the basic struc- 
ture of the magnetosphere. The strong corotating magnetic field of order 10 8 — 10 12 G 
dominates the physics in the surrounding region. The rotating magnetic field produces 
forces sufficient to rip electrons from the stellar surface with high energy. The high-energy 
electrons emit photons which in turn quickly decay via pair production. This process re- 
peats until the energies fall sufficiently to halt pair production. The consequent corotating 
charge distribution cancels the v x B forces, except perhaps at certain gaps, and the intense 
magnetic field constrains the electrons to move along the field lines. The end result is a 
one-dimensional relativistic electron-positron wind flowing out along the corotating mag- 
netic field lines. Those particles on the field lines which extend beyond the light cylinder 
can escape to infinity. 

Connecting this theoretical picture with observation has proven to be difficult; we still do 
not understand the physical mechanism which generates pulsar radio emissions. Observations 
of pulsars reveal that the individual pulses from a pulsar vary greatly from pulse to pulse 
with no apparent pattern. Averaging several hundred pulses, however, yields a stable profile 
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characteristic of a pulsar, yet these characteristic profiles still differ widely from pulsar to 
pulsar (Lyne & Graham-Smith 1998). Rankin (1983) has proposed that pulsars exhibit 
two distinct types of emissions, core and conal, suggesting that more than one emission 
mechanism may be at work. 

Despite the large variety of integrated pulse profiles, the profiles commonly exhibit a high 
degree of linear polarization. The polarization position angle often sweeps smoothly through 
a range of angles of up to 180°(Lyne & Graham-Smith 1998). In their rotating-vector 
model, Radhakrishnan & Cooke (1969) hypothesized that the radiation was polarized along 
or orthogonal to the magnetic field lines; the position angle therefore reflects the geometry 
of a dipolar magnetic field. 

While some pulsars exhibit an abrupt jump between orthogonal modes of polarization, 
their position angle distributions indicate that the emissions contain two orthogonal compo- 
nents, each of which appear consistent with the rotating-vector model (Manchester 1975; 
Backer, Rankin, & Campbell 1976). Fluctuations in the flux densities of the two modes 
would be observed as a sudden jump in position angle (McKinnon & Stinebring 1998). 

The pulse also typically broadens with decreasing frequency, a behavior usually at- 
tributed to radius-to-frequency mapping: lower-frequency emissions are thought to originate 
farther away from the stellar surface than higher-frequency emissions (see Cordes (1978) and 
references therein). Because the magnetic field lines diverge, the higher altitude translates 
to an increase in the size of the emission region which in turn translates to a wider observed 
pulse at lower frequency. 

Barnard & Arons (1986) investigated the propagation of electromagnetic waves in 
the relativistic plasma permeating the pulsar magnetosphere and found that one mode, 
the subluminal O-mode, follows paths along the magnetic field lines. A possible source of 
this radiation is the pair-production front, where the primary particles decay into those 
which make up the outflowing plasma (Lerche 1970; Melrose 1996). Lyutikov (1998) has 
also suggested that cyclotron instabilities in the plasma will produce subluminal O-mode 
radiation in the magnetosphere of a pulsar. 

Observations of rays originating in this mode would reflect the geometry of the dipolar 
field (Gallant 1996), which could account for some observed pulse shapes and position angle 
sweeps. In this paper, we consider a model that assumes a neutron star with a corotating 
dipolar magnetic field and a source of radio-frequency radiation near the stellar surface. By 
determining the path of rays from the stellar surface to the observer, we attempt to reproduce 
the pulse intensity and polarization profiles of a pulsar, specifically the Vela pulsar. 
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2. Wave Propagation in the Pulsar Magnetosphere 

In a relativistic plasma and a strong magnetic field B, electromagnetic waves propagate 
in two independent, linearly polarized modes. The electric field of an ordinary mode or 
O-mode wave lies in the (k, B)-plane, where k is the wave vector of the wave; the electric 
field of the extraordinary mode or X-mode points along the normal to this plane. Because 
the charges can move only along the field lines, the plasma has no effect on the X-mode, so 
this mode propagates as if in a vacuum. 

In the limit of an infinitely strong magnetic field, the dispersion relation for an O-mode 
wave of frequency uj propagating in a cold plasma with plasma frequency uj v and consisting 
of electrons and positrons moving with speed f3 p is given by 

where 7 P is the usual relativistic factor and k±_ and k u are respectively the components of the 
wave vector normal to and along the magnetic field (Arons & Barnard 1986). Figure 1 shows 
a plot of the dispersion relation in the rest frame of the plasma for various values of k' ± /u' p . 
Here the primes denote quantities measured in the rest frame of the plasma. For each value 
of k' ± /u' p , the figure shows two curves; the upper branch corresponds to the superluminal or 
fast mode while the lower branch corresponds to the subluminal or slow mode. 

Each ray begins its journey as a slow-mode wave and either transitions to the fast 
mode and escapes the magnetosphere or succumbs to Landau damping and dissipates in the 
plasma. We assume the transition occurs via some mechanism, when it reaches a point such 
that the distance between a pair of curves is minimized, or in other words, when the two 
types of waves are most alike. We approximate this criterion by requiring the transition to 
occur when k'Juj' p = 1. 

Several points should be made here. First, the criterion which determines where the ray 
escapes from the magnetosphere will obviously affect the observed pulse shape. The purpose 
of this paper, however, is not to explore the physics of the mode conversion but rather to 
examine how refractive effects coupled with the geometry of the dipole field might give rise 
to features of the pulse profile. Second, the transition criterion we chose corresponds to mode 
conversion in the weak-turbulence limit (Lyutikov 1998). Turbulence, in effect, smears out 
the dispersion curves, and the transition occurs via the region in which the smeared curves 
overlap. Third, the conversion rate may be so low in absolute terms so that it could not 
account for the observed radio emission. Clearly, the physics of the mode conversion is 
critical for a complete understanding pulsar radio emission and has been investigated (see, 
for example, Beskin, Gurevich & Istomin (1988)). However, such a discussion is beyond the 
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intended scope of this paper. 

We must also consider that the subluminal mode may interact with the electrons and 
positrons and transfer its energy to the plasma. The mode conversion can not occur if 
Landau damping dissipates the wave before it can reach its would-be transition point, but 
Arons & Barnard (1986) found that for canonical pulsar parameters, Landau damping would 
not be significant if the mode conversion occurred at radii less than 1000 km. We find that 
the transition points lie at altitudes significantly less than 1000 km and therefore neglect 
Landau damping. 

In terms of u' p 2 = uJ p /^/ p and uj 1 = 7 P (c^ — f3 p k u ) , where again the primes denote quantities 
measured in the rest frame of the plasma, the dispersion relation becomes 

(* 2 -*?)(i —I •■• ■■> (2) 




It follows that for the subluminal mode because u/ < uj' p , we must have |u;| < \k u \. In the 
corotating frame, the plasma is essentially moving at the speed of light and outruns the wave; 
boosting to the rest frame of the plasma therefore yields a negative frequency. Consequently, 
the relevant branch of the dispersion relation is given by 



« = ~\ V — (3) 

To determine how the wave evolves, we must know how k' ± and ui' p change as the wave travels 
away from the stellar surface. 

The slow O-mode wave packets obey the relations (Barnard & Arons 1986): 
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where r and 9 are the spherical coordinates of the wave packet in the corotating frame, with 
the dipole axis pointing along the illustrated in Figure 2. Here, the subscript 

denotes values at the stellar surface. The first equation simply says that the wave packet 
follows the magnetic field lines. The second equation describes how the direction of the wave 
vector relative to the field line evolves as the wave packet moves down the field line. As the 
wave moves farther away from the neutron star, the orientation of the wave vector depends 
less on its initial orientation and increasingly on its displacement from the star. 
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The strong magnetic field constrains the electrons and positrons to move along the field 
lines; hence, their density is proportional to the magnetic field strength which, for a dipolar 
magnetic field, varies as r~ 3 . In the corotating frame then, the plasma frequency is given by 

2 4vriV e 2 /r \3 

Up\ r ) = 2 — 5 

y m V r / 

where N is the plasma density at r = r , e is the elementary charge, and m is the mass of 
the electron. The factor of two reflects the fact that the plasma consists of positrons and 
electrons. Transforming to the comoving frame introduces a factor of 1/7 P because of length 
contraction, giving 

8nN e 2 

where uj' Pq = uj' p (r ). In reality, the magnetosphere may contain regions where uj p does not 
scale in this manner, and refraction may occur into these regions. Our model neglects this 
possibility. 

Knowledge of both k' ± and uo' p at each point on the trajectory allows us to determine 
which dispersion relation curve the wave satisfies but not the specific values of uj' and k' u . To 
find these quantities, we note that the wave also satisfies the Lorentz transformation 
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where u is constant in the geometric-optics limit. Plotted on the same axes as the dispersion 
relation, this relation corresponds to a line with slope — f3 p and intercept u/ r ) p oo l p . The 
intersection of this line with the dispersion relation curve gives us the values of uj' and k' u for 
the slow-mode wave at a given radius. Figure 3 shows four pairs of curves, corresponding 
to distances of ten, twenty, thiry, and forty stellar radii, for a ray following a field line. The 
dashed curve shows how u/ and k' n evolve as the ray travels away from the stellar surface. 

To find the transition point, we expand equation (3) in powers of k' ± /u' p with k'Ju' = 1 
to find 

uj' 1 k' ± 3 uj ( r \ 2 

— » -1 + = -1 + -—e - s 

u p 2^/ 16 uj' po \r J 

Similarly, from equation (7), we have 

w 1 „ uj uj ( r \ 3 ^ 2 . . 

— = -P+ — 7 = ~P + — r - 9 

IpUp lpUJ' po \r J 

Setting these two expressions equal to each other gives 

lpU' P0 VoJ 16 uj' P{) \r J 
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which may be solved numerically for r/r . 

After the mode conversion, a ray escapes the magnetosphere as a vacuum electro- 
magnetic wave. The angle between the wave vector and the field line is approximately 
n_i_ (3/8)0, and the angle between the field line and the dipole axis is (3/2)9. Therefore, 
after the transition occurs, the wave travels along a line that makes an angle £ = (9/8)0/ 
with the axis of the dipolar field, where 0/ is the value of at the mode transition. 

We simplistically assume in our model that the direction of the wave vector and the 
polarization remain unchanged by the mode conversion. If the two modes couple due to 
turbulence within the plasma, one would not expect such a conversion mechanism to preserve 
k and the polarization. Clearly, even in this simple model, the direction of the wave vector 
must change during the transition; because the phase velocities of the two modes are different, 
the constraints on uo and k± require that the two modes must have different k„s. The exact 
manner in which the slow mode maps to the fast mode depends on details of the conversion 
process, a discussion of which is beyond the scope of this paper. 

3. Ray Trajectories 

To analyze the path of a wave in the magnetosphere, it is convenient to focus on two 
reference frames: the observer's frame and the corotating frame. Figure 4 illustrates the 
relationship between cartesian coordinates in the two frames. The set denoted by capital 
letters corresponds to the observer's frame. The Z-axis points along the rotation axis of the 
pulsar, and the line of sight lies in the XZ-plane. The second set, denoted by lower-case 
variables, lives in the corotating frame. It is related to the first set by a rotation by i about 
the F-axis followed by a rotation by 5 about the Z-axis, where the angles S and i describe 
the orientation of the pulsar. In this system, the magnetic axis coincides with z-axis, and 
the rotation axis lies in the xz-plane. Finally, because we are working with a dipolar field, 
the corresponding set of spherical coordinates (r, 0, 0) in the corotating frame is useful. 

The transition from slow to fast mode may occur at a distance which is a significant 
fraction of the light cylinder radius; consequently, we must take into account the effects of 
the boost from the comoving frame to the observer's frame on the wave vector of the ray. 
The boost changes both the frequency and the direction of the ray. The frequency shift is 
small, about 1%, and is relatively unimportant. The change in direction, on the other hand, 
is more significant since the final direction of the ray in the observer's frame determines if 
and when the ray will reach the observer. 

Each point on the transition surface appears to emit a beam of light with a fixed direction 
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in the corotating frame. In the observer's frame, the beams rotate with the star. If a beam 
makes an angle a with the rotation axis, it will sweep across the line of sight when the star 
has rotated to an appropriate orientation, and the beam will reach the observer. 

Figure 5 illustrates the paths taken by several of these rays as seen from the corotating 
frame. Up to the transition point, the rays follow the magnetic field lines; after the transition, 
they travel in the directions given by their wave vectors. In the observer's frame, the rays 
will travel in straight lines after the transition. In the corotating frame, they will follow 
spiral trajectories. 



4. Constructing the Pulse Profile 

To construct the pulse profile, we start with a uniform distribution of rays on the stellar 
surface centered on the magnetic pole, trace their paths, and select only the ones which will 
reach the observer. The pulse period is divided into some number of equally-sized bins, and 
the time of arrival of a ray determines into which bin it goes. Averaging the position angles 
of the rays contributing to a bin yields the final position angle for that bin. 

Each ray carries a Gaussian weighting, with the distribution centered on 9 = 0. The 
width oq of the Gaussian corresponds physically to the size of the beam rising from the 
stellar surface. Appealing to the central limit theorem, we used a Gaussian distribution so 
that the resulting profile is not that of a single pulse but that of many averaged together. 
In other words, the model produces the integrated pulse profile. 

For a particular orientation of the neutron star, namely 5 = 0, we calculate for each ray 
its wave vector K in the frame of the observer. The time of arrival for a ray that will reach 
the observer is then given by 

tarr Q c W 

where 4>k is the azimuthal angle of the wave vector, Q is the rotational speed of the pulsar, 
R is the point where the ray breaks free of the magnetic field line, and K is the unit vector 
in the direction of the wave vector. The first term accounts for the fact that the star must 
rotate to the proper orientation for the ray to align with the line of sight. The second term 
arises from the displacement of the transition point from the star, as illustrated in Figure 6. 
The travel time for the light ray from a pulsar at a distance D is approximately t = D/c; 
the actual travel time differs from this value because the wave decouples from the plasma at 
a displacement R from the center of the star. 



As mentioned earlier, the rotating- vector model plausibly accounts for the characteristic 



- 8- 



sweep in position angle under the assumption that the polarization reflects the geometry of 
the dipolar field. The refraction model of wave propagation naturally explains the alignment 
of the electric field with the field lines because, for both the O-mode and X-mode, the 
direction of the electric field is determined by the directions of the magnetic field and the 
wave vector k. 

For a dipolar field, the rotating-vector model predicts that the position angle ip will 
vary according to 

, , rx sin a sin S _ 

tan V (5) = : 12 

sin l cos a — cos i sin a cos S 

where S = £lt. This relation, however, neglects relativistic and time-of-arrival effects. The 

transformation from the corotating frame to the observer's frame changes the direction of E. 

Moreover, because of the second contribution to the time-of-arrival delay, a ray corresponding 

to an orientation Si can arrive before, at the same time as, or after a ray corresponding to 

S2 > Si. These effects may help explain the deviations observed by Krishnamohan & Downs 

(1983) from the predictions of the basic rotating-vector model. 

We obtain the polarization of a ray in the observer frame by calculating the electromag- 
netic field tensor F^ u for the wave in the corotating frame and then transforming it to the 
observer frame. To find the electric and magnetic field components in the corotating frame, 
we define a set of cartesian coordinate axes so that the z-axis points along the wave vector 
and the x-axis lies in the (k, B)-plane. The electric field then points along rotated rr-axis, 
and the magnetic field, along the y-axis. 

We have assumed that the electric field of a ray lies in the plane of the dipole field line. 
Cheng & Ruderman (1979) found that a narrow cone of rays emitted about a field line with 
initially no average polarization will eventually become nearly completely polarized and that 
the electric fields will lie in the plane of the dipole field line, if the cone and the direction of the 
magnetic field become sufficiently misaligned. For radiation due to longitudinal acceleration 
parallel to B, the opening half angle of the cone is of order 7" 1 0.01 rad. The angle 
between the cone and the field line at the transition point is of order Of m 0.1 rad; hence, 
our assumption seems reasonable. 

5. Results From Applying The Model To The Vela Pulsar 

Our model has eight parameters. Six of these describe the pulsar: i?* is its radial size; 
T*, its rotation period; L*, its spin-down luminosity; 1, the angle between the magnetic field 
axis and the rotation axis; j p , the bulk relativistic factor of the electron-positron wind; and 
ae, the width of the beam. The two remaining parameters a, the angle between the rotation 
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axis and the line of sight, and ou, the observing frequency, describe the observer. 

All of the parameters are set by observations, theory, or the observer except for a, t, 
00, and 7 P . However, observations of polarization often determine the quantity a — i, and 
sometimes determine i as well. Thus, our model contains, in effect, two or three parameters 
that can be set to change the pulse shape. 

The degree to which refractive effects dominate the propagation of the electromagnetic 
waves depends mainly upon the plasma density on the open field lines. We calculated this 
density by assuming that the spin-down luminosity of the pulsar is carried away by the 
electrons and positrons on the open field lines; thus, the ratio L*/7 P sets the density of the 
electron-positron wind at a given radius. We chose to apply the model to the case of the 
Vela pulsar because its relatively high spin-down luminosity = 6.9 1 8 x 10 36 erg s _1 and 
correspondingly dense plasma should emphasize the refractive effects. 

Vela's short 89.30 ms rotation period translates to a light cylinder radius of 4236 km. 
We found that the transitions from slow-mode to fast-mode waves occur within 10% of this 
distance. Relativistic effects, small as they are, could still have a significant effect on the 
pulse shape. 

We set i to 90° and to 10 km and chose the values for the remaining parameters 
7 P = 100, oq = 0.020 rad, and a = 94.8° so that the model would produce a reasonable pulse 
shape for an observing frequency of 2.3 GHz. We did not perform an extensive search of the 
parameter space to study the effects on the pulse profile. 

Figure 5 depicts the paths of several rays to the observer as seen in the corotating 
frame. In this frame, observers follow paths which encircles the magnetic pole and see a 
pulse when they cuts across the "spray" of light. From this figure, we also see that observers 
would measure the size of the pulsar, i.e. the size of the transition surface, to be around 
100 km at 2.3 GHz, a result which agrees reasonably well with recent size measurements of 
Vela. Macquart et al. (1999) obtained an upper limit of 50 km from 660 MHz scintillation 
data and suggested that the size of the emission region increases with frequency; Gwinn et 
al. (2000) measured the size of Vela's emission region to be 440 km; and Krishnamohan & 
Downs (1983) found a spread in altitude of about 400 km. 

Figure 7 shows the intensity profiles obtained from the model for two frequencies and 
the associated transition surfaces. The dotted curves correspond to 1.4 GHz, and the solid 
curve, to 2.3 GHz. The larger transition surface at 1.4 GHz results in a wider pulse. A plot 
of the position angle shows a smooth sweep, somewhat obscured by the break in the curve 
caused by the coarse graining of the bins. At the start of the pulse, the shape of the curve 
deviates from the S-shaped curve predicted by the simple rotating- vector model because rays 
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from many longitudes contribute to the leading edge of the pulse. 

Figure 8 shows the observed integrated pulse profile for the Vela pulsar at 1.4 GHz. 
The average position angle sweeps through roughly the same range for both profiles. More 
interesting, the predicted and observed profiles both appear to consist of the overlap of a 
strong subpulse followed by a weaker subpulse, although the model did not reproduce the 
relative strengths accurately. In our model, the presence of two subpulses arises from a 
combination of three effects. 

Every ray reaching the observer emanates from a point on the stellar surface; therefore, 
we can map a pulse to a curve or a region, if rays from multiple points reach the observer at 
the same time, on the stellar surface. Moving along this curve corresponds to advancing in 
time on the profile. Typically, this path starts at some point (#o> 0o) — (©o>^o) an d while 
(p monotonically decreases, 9 first decreases, reaches a minimum, and then increases. 

Because of the cylindrical symmetry around the magnetic pole, each circle of latitude 
on the stellar surface has an associated circle of transition points in the magnetosphere. The 
density of rays reaching the observer is higher when the line of sight grazes one of these 
circles than when it cuts across the circle at some angle. Hence, the density of rays reaches 
a peak when do hits its minimum and the path on the stellar surface is tangent to a circle of 
latitude. 

Each ray represents an area element dVL$ ~ 6od6od<po, so its weight contains a factor 
of #o as d9o and d<fio are fixed in our model. The weight assigned to a ray also has the 
Gaussian factor mentioned earlier, representing the strength of an average beam. The beam 
is strongest at the magnetic pole and weakens away from the pole with a characteristic 
distance specified by a e . 

This combination of weighting factors can easily lead to the appearance of two subpulses. 
At first, when 6 is large, the Gaussian factor suppresses the weighting. As 9 decreases, the 
weighting at first steadily increases. At some point, the area factor 6q will begin to dominate, 
and the weight will decrease, resulting in the first subpulse. After 9 reaches its minimum, 
the reverse sequence occurs, generating the second subpulse. 

Neither the higher density of rays at smaller values of #0 nor the weighting assigned to 
the rays can result in the observed pulse with its asymmetry between the two subpulses. 
Neither effect nor their combination breaks the cylindrical symmetry. If the density effect 
enhances one weighting maximum, it will enhance the other one equally. With only these 
two effects, the observer would always see symmetrical pulse. 

The final effect is the time-of- arrival effect, which does break the symmetry. In the 
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corotating frame, the cylindrical symmetry about the magnetic pole results in a transition 
surface with a reflection symmetry through the x-z plane. The Lorentz boost to the observer's 
frame, however, breaks the symmetry by rotating all of the wave vectors in the direction of 
the boost and, because this rotation affects the time-of-arrival, skews the pulse to the left in 
the profile. The other contribution to the time-of-arrival delay further skews the pulse to the 
left; rays from further out on the transition surface take less time to travel to the observer 
and arrive earlier in the pulse. Also because their transition points are moving faster, they 
undergo a larger boost, resulting in a larger rotation and causing them to arrive even earlier 
in the pulse. 

Figure 9 illustrates how these effects combine to give the observed pulse. The top plot 
shows that rays starting near the magnetic pole arrive early in the pulse. The bottom plot 
shows the weighting with two maxima, but these small variations can not alone account for 
the large difference in strength between the two subpulses. The superimposed histogram 
shows that many rays amass near the beginning of the profile. The first subpulse results 
then from a combination of the first maximum in the weighting and the large number of rays 
reaching the observer at the same time. The second subpulse arises as the rays arrive at a 
constant rate but with decreasing weights. 



6. Summary 

A simple ray-tracing model that assumes a neutron star with a corotating dipolar mag- 
netic field yields results which appear qualitatively consistent with observations. The re- 
sulting pulse profile exhibits the general characteristics of the observed profiles; the size of 
the emission region agrees with recent measurements; and the model naturally incorporates 
the rotating-vector model if one assumes the mode conversion process does not obliterate 
polarization information. 

The basic nature of the model, moreover, leaves room for refinement. We have neglected 
the physics in the transition from the slow mode to the fast mode, and an improved model 
may involve a more realistic momentum distribution for the plasma. Finally, adjusting 
the weight assigned to each ray by incorporating a more detailed model of the radiation 
distribution at the surface may yield information about the shape of individual pulses and 
their variations. 
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A. Coordinate Transformations Prom the Corotating Frame to the Observer's 

Frame 

After solving equation (10) for r/r , we can find Of using equation (4a). Due to the 
cylindrical symmetry of the dipole field, transitions occur for all values of 0. For a given 
value of 0, the transition occurs at the point 

r = r(sin6>/ cos0, sin6>/ sin0, cos Of), (Al) 

and the ray has a wave vector 

k = u;(sin£cos 0, sin £ sin 0, cos£) (A2) 

where £ = (9/8)0/. The electric field e and magnetic field b are given by 

e = (cos £ cos 0, cos £ sin 0, — sin£) (A3a) 
b = (— sin 0, cos 0,0). (A3b) 

We are working in units where c = 1. The vectors are normalized as we are not con- 
cerned with their magnitudes but only with their directions so that we can determine the 
polarization of the ray as seen by the observer. 

The coordinates in the corotating frame and the observer's frame are related via two 
rotations so that 

R = (X,Y,Z)=R z (8)Ry(i)t (A4a) 
K = (K x ,K Y ,K z )=R z (5)R Y (i)k (A4b) 
E = (E x ,E Y ,E z )=R z (5)R Y (L)e (A4c) 
B = (B x ,B Y ,B z )=R z (5)R Y (L)b (A4d) 

where R z and R Y are the rotation matrices for rotations by the specified angles about the 
Z-axis and F-axis respectively. 



B. Lorentz Transformations 

The vector K does not give the wave vector of the ray as seen by the observer; rather, 
it is the representation of the wave vector as measured in the corotating frame expressed in 
the coordinates of the observer's frame. To find the wave vector as seen by the observer, 
we must apply a Lorentz transformation because the transition point moves with a velocity 
(3 = (ftZ) x R/c. 
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The components of K parallel and perpendicular to j3 are given by 

K„ = (K-P)P (Bla) 
K ± = K-K, (Bib) 

where (3 denotes the unit vector in the direction of (3. According to the observer, the ray 
has a wave vector K = K n + Kj_ where 

K = <y(K H + l3u)P (B2a) 
K ± = K ± . (B2b) 

and 7 = l/y/l — (3 2 is the usual relativistic factor. 

We must likewise transform the electric and magnetic fields. The electric and magnetic 
fields that the observer will see are given by 

E = 7 (E-/3xB)-^-(/3.E)/3 (B3a) 

7 + 1 

B = 7 (B + /3xE)-^-(/3.B)/3. (B3b) 

7 + 1 

The position angle ip is then given by 

■0 = arctan I = 1 (B4) 



^E 2 - E 
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Fig. 1. — Plot of the cold plasma O-mode dispersion relation in the rest frame of the plasma 
for various values of k' ± /u' p , as indicated by the numbers labeling each curve. The upper 
branch corresponds to the superluminal or fast mode; the lower branch, to the subluminal 
or slow mode. 




Fig. 2. — As shown by the halftone line, the slow-mode wave leaves the stellar surface at 
(R*,9o) and follows a magnetic field line. At (rf,0f), it transitions to the fast mode and 
travels along kf. 
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Fig. 3. — A line and a dispersion relation curve are associated with each point in the 
trajectory of an O-mode wave packet. The intersection of the two yields the waves values of 
uj' and k' u at that point. The dashed line shows how uj' and k' n evolve as the wave propagates 
away from the pulsar. 




Fig. 4. — The axis of rotation points along the Z-axis. The line of sight lies in the X, Z-plane 
and makes an angle a with the rotation axis, and the direction of the magnetic axis is given 
by the polar angle i and the azimuthal angle 5 = Qt. 
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Fig. 5. — These figures show the paths of several rays, as seen in the corotating frame, 
projected onto the xz- and yz-planes. 




Fig. 6. — The ray leaving from point P travels a distance R • k less than a ray from the 
origin. Hence, it will arrive at the observer a time 5t — R • k/c earlier. 
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Fig. 7. — This figure shows model-generated pulse profiles at 1.4 GHz (dotted) and 2.3 GHz 
(solid) and the corresponding transition surfaces. The pulse period consists of one hundred 
bins. 




Fig. 8. — The observed integrated pulse profile of the Vela pulsar at 1.4 GHz. We thank S. 
Sallmen and D. C. Backer for providing these data. 
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Fig. 9. — Rays from near the magnetic pole arrive early in the pulse, enhancing the first 
subpulse. 



